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Constraint on intergalactic dust from thermal history of 
intergalactic medium 
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ABSTRACT 

This Letter investigates the amount of dust in the intergalactic medium (IGM). 
The dust photoelectric heating can be the most efficient heating mechanism in the 
IGM where the density is very small and there are a lot of hard ultraviolet photons. 
Comparing the observational thermal history of IGM with a theoretical one taking 
into account the dust photoelectric heating, we can put an upper limit on the dust- 
to-gas ratio, T>, in the IGM. Since the rate of the dust photoelectric heating depends 
on the size of dust, we find the following results: If the grain size is > 100 A, V at 
z ~ 3 is < 1/100 Galactic value corresponding to fiJuirb ~ 10~ 5 - O n the other hand, 

if the grain size is as small as ~ 10 A, T> is < 1/1000 Galactic value corresponding to 
QlGM < in-6 
"dust ~ 1U ■ 
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1 INTRODUCTION 

There are metals in the intergalactic medium (I GM) even in 
low density regions such as Lyman a forest (e.g.. lCowie et alJ 
Il995l: ITelfer et aljliooj) . Since metal and dust are relating 
each other, it is sure that dust grains also exist in the IGM. 
However, the amount of the intergalactic (IG) dust is still 
quite uncertain although it seems to be not so abundant. We 
try to put a new constraint on the amount of the IG dust 
by using the IGM thermal history s uggested by the rec ent 
observations of Lyman a forest (e.g.. ISchave et alJl2000l) . 

After the early attempts to estimate the IG extinc - 
tion (lEigensonl Il949t iHumason. Mavall. fc Sandagdfl95d) . 
ICrane fc Hoffman! J1973TT have obtained an upper limit on 
the amount of the IG dust in terms of the density parameter 
as fiJiu^ = 0) < 10 -4 , by comparing the observed spectral 
energy distribution of distant (but z < 0.5) giant elliptical 
galaxies with an average spectrum of similar nearby galaxies 
(see also iNickerson fc Partridgeill97lt lTakaselll972l) . 

Observations of the redshift evolution of the QSOs spec- 
tral slope provide us with another approach to t he IG dust 
JWrighdll98ll : ICheng. Gaskell. fc Koratkarlll99ll) . Although 
Cheng et al. (1991) found no evidence for an appreciable 
IG reddening up to z ~ 2-3, this may not be so strict con- 
straint because it is difficult to measure the UV slope of 
QSOs enough precisely owing to the broad and complex iron 
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band superposing on the continua. Recently, we obtain fur- 
ther constraint on the IG dust from observations of high- z 
supernovae (SNe) llRiess et aJl99SUPerlmutter et al.ll999l ). 
An upper limit of the colour excess by the IG reddening is 
(E B -v)z~a.s - { E B -v)z~am < 0.03 mag. 

By the way, lAguirre) l)l999l) suggests that the IG dust 
can be "gray" by a selection rule in the transfer of grains 
from the host galaxies to the IGM. This occurs when small 
(< O.lum) grains are destroyed selectively by the thermal 
sputtering in the hot gas halo of the host galaxies before 
they reach the IGM. Then, the extinction property becomes 
"gray" since the relatively large dust survives. Importantly, 
such "gray" IG dust can avoid the detection by the IG red- 
dening survey like above. 

Even if the IG dust is really "gray" , its evidence should 
be imprinted in the cosmic microwave background (CMB) 
and infrared background because the dust emits thermal 
radiation in the wav e-band from the far-infrared to sub- 
millimetre (submm) l|Rowan- Robinson. Negrooonte. fc Silkl 
Il979t IWrighti Il98ll) . According to the current status, 
CO BE data provides us with a rough upper limit on 
the IG dust llLoeb fc HaimarJ Il997t iFerrara et all Il999l : 
lAguirre fc HaimarJ 12000 ) . On the other hand, the submm 
background radiation may give a more strict constraint on 
the IG dust. The emission from the IG dust expected from 
the modern cosmic st ar formation h istory with the grain 
transfer mechanism of lAguirrel (Il99flt) contributes to a sub- 
stantial fraction (> 75% ) of the measured back ground radi- 
ation in the submm jAguirre fc HaimarJ Eoool) . It is inter- 
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esting to develop the approach of lAguirre fc HaimarJ J200(J) 
to match a strong constraint by SCUBA whose result of the 
numb er count accou nts ~ 90% of the submm background 
light jCalzettill200ll and references therein) . 

In any current status, the amount of the IG dust is still 
highly uncertain. Hence, it is worth trying to find a new 
constraint on the amount of the IG dust. In this Letter, it 
is demonstrated that we can obtain such a new constraint 
of the IG dust by comparing the observational IGM tem- 
perature with a theoretical one if we take into account the 
dust photoelectric heating in the theoretical model. This is 
possible bec ause the dust heating is an effici ent mechanism 
in the IGM <Nath. Sethi, fc Shchekinovll999l) . Interestingly, 
the effect of the dust photoelectric heating depends on not 
only the amount but also the size of dust. Therefore, we can 
put a limit on the amount of the IG dust as a function of 
the typical size of the IG grains. 



2 PHOTOELECTRIC HEATING BY GRAINS 

To assess photoelectric effect, we must specify the charge on 
grain, Zd (in the electron charge unit), which is given by 
JSpitzedll94ll) : 



^ Ri + Rp 



0. 



(1) 



where we have assumed a charge equilibrium because a typ- 
ical charging time-scale is very short (~ 10[a/0.1um] _1 yr 
under the UV background radiation dominated by QSO light 
given in §3). In the equation, Ri is the collisional charging 
rate by i-th charged particle, and R pe is the photoelectric 
charging rate. We assume spherical grains for simplicity. 
The collisional charging rate is 



Ri = na ZiSiUi 



i(vi, Z d , Zi)vif(vi)dvi 



(2) 



where Zi is the charge in the electron charge unit, Si is 
the sticking coefficient, m is the number density, a is the 
grain radius, Vi is the velocity, w m i n is the minimum velocity 
required to collide with a grain, en is the dimensionless col- 
lisional cross section depending on both charges and Vi, and 
f(vi) is the velocity distribution function. We simply assume 
Si is always unity. In the above equation, we neglected the 
effect of the secondary emission. 

We can neglect the effect of the "image potential" 
JPraine fc Sutinlll987l) on the collisional cross section be- 
cause the obtained charges are enough large. If we as- 
sume the Maxwellian velocity distribution for the parti- 
cle, the integral in the r.h.s. of equation @ is reduced to 
(8k B T/Tvmi) 1/2 g(x) and g(x) = 1 - x for Z d Z t ^ or 
g(x) = exp(— x) for Z^Zi > 0, where fee is the Boltzmann's 
constant, T is the gas temperature, m; is the charged par- 
ticle's mass, and x = e 2 Z^Zij 'aksT . 

The photoelectric charging rate is given by 



Rn 



Q(a, u)Y(a, v, Zd) ^^J dv , 



(3) 



where Q is the absorption coefficient of grains, Y is the pho- 
toelectric yield, J„ is the intensity (averaged for solid angle) 
of the incident radiation at a frequency v, h is the Plank con- 
stant, and i^max is the maximum frequency of the incident ra- 



diation. The minimum frequency, Vmin, is the threshold pho- 
ton frequency required for an electron to escape from a grain. 
That is, hfmin is equal to the ionization potential, IP, of a 
grain: IP = W + (Zd + l/2)e 2 /a, where W is the work func- 
tion. Here we neglect quantum effe cts on IP because it is 
very s mall for grains with a > 10 A iWeingartner fc Draind 
l2001al. hereafter WD01). 

According to WD01, we adopt W — 4.4 eV for graphites 
and W = 8.0 eV for silicates. For Q, we adopt the val- 
ues of "graphi te" and "smoo t hed UV astronomical s i licate " 
calculated bvlDraine fc Led (Il984l) ; lLaor fc Draind Jl993>) : 
IWeingartner fc DraindT^OOll ). The value of Y is estimated, 
based on the way constructed by WD01, which includes ap- 
proximately the effect of the energy distribution of photo- 
electrons and the geometrical enhancement for small grains, 
and reproduces recent results of the laboratory experiments 
even for very small grains. 

The both equilibrium charges obtained for graphites 
and silicates are similar when their radii are equal. This 
is consistent with the result of Nath et al. (1999). In the 
following discussion, we consider mainly silicate case. 

The photoelectric heating rate by a grain with radius a 
is expressed by 



7(a) 



E pe (a, v, Z d )QY ^ Jv dv , 
hv 



(4) 



and the cooling rate of the electron capture (recombination) 
by the grain is A(a) = (3/2)A;bT|_R c |, where we have assumed 
the Maxwellian distribution for the gaseous electrons. We 
also have defined E pa , the mean kinetic energy of photo- 
electrons that are emitted from a grain with radius a and 
charge Z d when a photon with energy hv is absorbed. This 
is determined from the energy distribution function of the 
photoelectrons, f(E), which is assumed to be a parabolic 
function introduced by WD01, 1 and is consistent with a 
typical energy distribution of photoelectrons in the labora- 
tory experiments. Hence, we estimate 

= / Ef(E)dE = ^" (£m " 7^ m ' n) , (5) 
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3E 

mm 



where S max = hv — IP is the maximum energy of the photo- 
electrons and £mim = —e 2 (Zd + l)/a is the minimum energy 
of the photoelectrons appearing on the grain surface (how- 
ever falling back into the grain). For the parabolic function 
adopted here, a typical kinetic energy of the photoelectrons 
is about half of -E max . 

To estimate the total photoelectric heating rate, we 
specify a characteristic grain size in the IGM instead of the 
size distribution of the IG grains. This is because the size 
distribution in the IGM is quite unknown. Here we consider 
three cases for the typical grain size; 0.001 nm, 0.01 u.m, 
and 0.1 u.m. The smallest case of 0.001 u.m is based on the 
size of the primary grai ns produced in the ejecta of SNe II 
jTodini fc Ferraral 120 011. It is interesting to obtain an in- 
formation about the size of the IG dust as well as to put 
limits on the amount. From the analysis of the case of very 
small IG dus t, we can asses s the selection rule of the IG dust 
suggested by lAguirrd dl99 



Our f(E) is equal to f E (E) = f%(E)/y 2 in WD01. 
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The grain number density is given as n(a) = p^/m(a), 
where m(a) = (47r/3)pa 3 is the mass of a grain with the 
radius a, g (= 2gcm~ 3 ) is the grain material density, and 
pd is the volume grain mass density, which is given by pd = 
pmnribD, where p is the mean atomic weight, ran is the 
proton mass, rib is the baryon number density, and T> is the 
dust-to-gas mass ratio. We introduce a relative dust-to-gas 
ratio, £, defined as D — CPm w. where T>m w — 6 x 10~ 3 is 
the Galactic dust-to-gas ratio ( Spi tzerll978f) . for convenience 
below. Therefore, the total photoelectric heating rate and 
electron capture cooling rate per unit volume are r pc = 
"f(a)n(a) and A pc = X(a)n(a), respectively. 

The total heating rate is approximately proportional 
to a~ ' 5 under the situation considered in the next section. 
This is because the heating rate per a grain roughly has 
a dependence of a 2 5 , i.e., the grain cross section plus an 
additional a dependence caused by J E pB QYdv, whereas the 
grain number density is proportional to a~ 3 . Therefore, for 
a fixed dust mass, the photoelectric effect becomes more 
important as the grain size is smaller. The cooling rate by 
the electron capture is always about an order of magnitude 
smaller than the photoelectric heating rate. 



3 DUST AND THERMAL HISTORY IN IGM 

Suppose an ideal fluid element with the mean density of the 
IGM. The time evolution of its temperature is governed by 
fe.g.. lHui fc Gnedinlll997l) 



A) 



(6) 



dT rr T dX 2(r 

2_f/T -|- 

dt X dt SksXrih 

where H is the Hubble constant, rib is the cosmic mean 
number density of baryon, T and A are the total heating 
and cooling rates per unit volume, respectively, X is the 
number ratio of the total gaseous particles to the baryon 
particles, i.e., X = nj/nh, where n; is the number den- 
sity of the i-th gaseous species and we consider H I, H n, 
He I, He II, He III, and electron. Then, we solve equation © 
coupled with the non-equilibrium rate equations for these 
gaseous par ticles . For the rate coefficients, we use those 
compiled bv lCenl jl992T) . The adopted initial conditions are 
T = T C mb(1 + z), H i/H(total)=0.99, He i/He(total)=0.99, 
He n/He(total)=0.0099 at redshift z = 10, where T C mb = 
2.7 K is the temperature of the CMB radiation at z = 0. 
The result is insensitive to the choice of the initial condition 
if calculation is started at an epoch before H I reionization. 
We also adopt H = 70 km s" 1 Mpc" 1 , fi M = 0.3, fi A = 0.7, 
fib = 0.04, and He mass abundance Y = 0.24. 

To add to the photoelectric heating and the electron 
capture cooling by grains, we consider the following heat- 
ing/cooling mechanisms: For cooling, we use the recom- 
bination cooling, collisional ionization/excitation cooling, 
brcmsstrahlung cooling, and Compton cooling compiled by 
ICenl jl992f) . For heatin g rates, to take accou nt of the ra- 
diative transfer effect llAbel fc Haehneltlll999l) . we multiply 
usual photoionization heating rates by a correction factor of 
(1 + Crt/Oi where fi is the fractional abundance of H I, He I, 
and He II relative to total number of H and He. This factor 
mimics the situation that atoms are ionized by higher en- 
ergy photons when the optical depth is large (i.e., /, is nearly 
unity) relative to the optically thin limit. The coefficient Crt 



is a parameter determined by solving the cosmological ra- 
diative transfer. The metal line cooling is negligible because 
we consider only the case that metallicity is less than 1/100 
sola r value and temperature is less than a few 10 4 K (see 
e.g.. lSutherland fc DopItalll993l) . 

We need to give the UV background radiation and the 
cosmic reionization history to obtain the temperature evo- 
lution. We adopt, for UV background, J v = Ji\v~ a , where 
J%\ is the mean intensity at the Lyman limit of H I in unit 
of 10 -21 erg s _1 c m" 2 sr -1 Hz" 1 and its evolution is given 
by equation (3) in iKitavama et alJ feOOlft . We assume the 
spectral index to be a = 1 as a QSO dominated back- 
ground radiation. For the reionization history, we assume, 
for simp licity, a sudden re ionization of H I (and He i) at 
z = 6.0 (iBecker et alJl200 ll) 2 : there is no UV photon before 
this redshift. To m imic the He II reionization at z » 3.4 
jTheuns et al]l20 02*) , we set the maximum photon energy of 
the UV background to be 54.4 eV (He II Lyman limit) for 
z > 3.4 and to be 1.24 keV 3 for z < 3.4 (i.e., a sudden He II 
reionization). Other choice of a will be discussed later. 

We compare our theoretical thermal history at the mean 
density of th e IGM with the observational one obtained by 
ISchave et alJ |2000). They observe the Lya forest clouds 
with the column density of 10 13 ~ 15 cm" 2 (i.e., slightly over 
density regions), and then they convert the temperature at 
this density range estimated from the minimum b (Doppler) 
parameters into that at the mean density of the IGM by 
using the observed equation of state of the IGM. Thus, we 
can compare both thermal histories directly. 

In Figure 1, we show the IGM thermal history obtained 
for no dust cases and dusty cases of £ = 1/100 (silicate 
grains). For no dust cases (three dotted curves), the calcu- 
la ted thermal h i stories are consistent with the data points 
of ISchave eta!] feOOOft if C RT ^ 3.0-7.0. For dusty cases 
(solid [a = 0.1 vim], dash-dotted [a = 0.01 vim], and dashed 
curves [a — 0.001 Vim]), we observe, especially after the He II 
reionization, temperature enhancement by the dust photo- 
electric heating. In Table 1, we summarize the temperature 
enhancement factors at some redshifts for some sets of dust 
size and dust-to-gas ratio. All calculations do not include 
the evolution of the dust-to-gas ratio, for simplicity. 

We also present the time evolution of heating rates per 
a unit volume for the case of a = 0.01 vim and £ = 1/100 in 
Figure 2. The solid, dashed, dotted, and dash-dotted curves 
are the dust photoelectric, H I, He I, and He II photoion- 
ization heating rates, respectively. These heating rates show 
sudden change at the He II reionization epoch (z = 3.4). 
This is because the maximum photon energy increases from 
54.4 eV to 1.24 keV suddenly at that time. We note here 
that the heating rate by grains is saturated even when the 
maximum photon energy increases more than about 300 eV 
in the case of the spectral index of the incident radiation 
of a = 1. Thus, our results are robust for the choice of the 
maximum energy after the He II reionization as long as we 
adopt the energy larger than 300 eV. In addition, we show 



2 Although Ko gut et alJ 120031) find a higher-z H I reionization, 
our constraint of the IG dust dose not change because it is deter- 
mined mainly from the IGM temperature after the He II reion- 
ization. 

3 This is determined by the maximum frequency in the data of 
grain absorption coefficient. 
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Table 1. Temperature enhancement by IG dust. 



Figure 1. IGM Thermal history under sudden reionizations of 
H I (and He i) at z = 6.0 and He II at z = 3.4. Data points with 
errorbars are taken from Schaye et al. (2000). Three dotted curves 
are no dust cases; bottom: optically thin case, middle: Crt = 3.0 
case, top: Crt = 7.0 case, where Crt is the correction factor of 
the radiative transfer effect. The dashed, dash-dotted, and solid 
curves are cases of grain radius a = 0.001, 0.01, and 0.1 um, 
respectively. For these dusty cases, the dust-to-gas mass ratio is 
set to be 1/100 Galactic value, Crt = 3.0 is assumed, and the 
grain type is silicate. The middle dotted curve in z > 3.4 is almost 
superposed on the solid curve. 
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Figure 2. Time evolution of heating rate per a unit volume for 
the case of grain radius a = 0.01 um, dust-to-gas mass ratio of 
1/100 Galactic value, and radiative transfer correction Crt = 3.0. 
Solid curve is the photoelectric heating rate by grains. Dash, dot- 
ted, and dash-dotted curves are the photoionization heating rates 
by H I, He I, and He II, respectively. Two spikes at z = 6.0 and 3.4 
are produced by H I and He I reionization, and He II reionization, 
respectively. Two thin solid curves are the photoelectric heating 
rate of the cases of a = 0.001 um (top) and a = 0.1 um (bottom). 
Grain type is silicate. 



the cases of a = 0.1 and 0.001 \im (and ( = 1/100) by thin 
curves in Figure 2. Even if we adopt other sets of a and (, 
the other curves of the photoionizations do not change so 
significantly. 

Let us discuss each case in more detail. For the cases of 
a = 0.01-0.1 u.m, before the He II reionization, the photo- 
electric heating by grains is minor heating mechanism rela- 
tive to the H I photoionization heating. Thus, the effect of 
the photoelectric heating on the IGM temperature is very 
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Table 2. Amount of the IG dust at z > 2. 
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small, i.e., temperature enhancement is less than a few %. 
After the He II reionization, the photoelectric heating dom- 
inates other heating mechanisms (but by a factor of less 
than 3) when £ = 1/100. However, the temperature enhance- 
ment factor is still less than 1.5 (solid and dash-dotted curve 
in Figure 1). Thus, the thermal histories of these ca ses are 
consistent with that obtained bv lSchave et alj (|2000) when 
Crt = 3.0 is adopted. 

For the case of very small (~ 0.001 u.m) grains, the 
photoelectric heating is always dominant when £ > 1/100. 
Especially, this heating exceeds the He II photoionization 
heating by a factor of 10 after the He II reionization. Then 
the temperature enhancement factor becomes about 2 at 
2 = 2 (dashed curve in Figure 1). Thus, the thermal his- 
tory in this c ase m ay be inconsistent with the data from 
ISchave et aT] fcOOOl) . This means that if the typical grain 
size in the IGM is ~ 0.001 [im, the dust-to-gas mass ratio 
in the IGM should be much less than 1/100 Galactic value. 



4 DISCUSSION 

In Table 2, we summarize the obtained upper limits of the 
amount of the IG dust at z > 2. Let us compare these limits 
with the IG dust amount es timated from a possible cosmic 
star formation history (e.g., iMadau. Pozzetti. fc Dickinson! 
1998). Once we assume a cosmic star formation history, the 
cosmic metal amount can be estimated. Roughly, we find 
fimctai ~ a few xlO" 5 at z ~ 3 (e.g.. lAguirreUl999l) . If 
we postulate ~ 0.5 as the metal depletion to dust grains, 
which number is suitable for the Milky Way, we obtain 
f^dust ~ 10 -5 as the total cosmic dust amount at the red- 
shift. We show in Table 2 the upper limits in terms of the 
density parameter converted from the dust-to-gas ratios by 
fidurt 0) = (fMw(2b(z). For the lar ger grain models, we ob- 
tain an upper limit, 



7 x 10 which is comparable 
to fidust- This means that our upper limit may not deny the 
possibility that most of the dust grains produced in galaxies 
before z ~ 3 escape to the IGM as an extreme case. For very 
small (0.001 u.m) grain case, we obtain fijgirt i$ 7 x 1(r? - 
Thus, a small fraction (< 10%) of such very small grains 
can escape from the host galaxies to the IGM because of 



; /^dust < 0.1. The latter conclusion is consistent with 
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the suggestion bv lAguirrel <ll999l) . who proposes that the 
IGM is very difficult to be polluted by the very small grains. 

The obtained values of the heating rate by the dust 
photoelectric effect are consistent with those in Nath et 
al. (1999). It indicates that our result is robust against 
uncertainties of the adopted efficiency factor Q, photoelec- 
tric yield Y, and energy distribution function of the pho- 
toelectrons, etc., because these are somewhat different be- 
tween our model and Nath et al. (1999). If we consider only 
graphite grains, the photoelectric heating rate is found to be 
about 70% of that of only silicate grains, i.e., silicate grains 
can be more efficient heating source. The temperature in the 
graphite case is lower than that in the silicate case, while the 
decrement is < 10%. These results are also consistent with 
Nath et al. (1999). 

We shall comment some uncertainties of our analysis. 
First, the effect of the spectral index of a is examined. As 
shown in Zheng et al. (1997), for example, a softer spectrum 
than a = 1 is still compatible with observations. If a — 1.5, 
the heating rates by dust at z — 3 decrease by a factor 
of 0.25 (a = 0.1 /im)-0.5 (a — 0.001 fim) relatively to our 
case of a = 1. If a = 2 is used, the heating rates at z — 3 
decrease by a factor of 0.1 (0.1 /im)-0.25 (0.001 /im). In these 
calculations, £ = 1/100 and Crt = 3.0 are assumed. As the 
dust heating becomes less efficient, a more abundant IG dust 
is allowed. Then the upper limits for the IG dust amount 
increase by a factor of 2 (0.001 /im)-4 (0.1 /im) for a = 1.5 
and by a factor of 4 (0.001 /im)-10 (0.1 jim) for a = 2 
relative to the a = 1 case, because the dust heating rate is 
proportional to the dust-to-gas ratio linearly. To summarize, 
a softer spectral index constrains more loosely the amount 
of the IG dust by a factor. 

It might be possible to take very large a. For example, 
we consider a background radiation field dominated by only 
galaxies (e.g., a = 5). In this case, the dust heating rates 
at 2 = 3 become 1/100 (0.1 /im)-l/10 (0.001 /im) of those 
in a = 1. However, the thermal histories with a = 5 never 
reproduce a steep rise of the observed temperatures at the 
He II reionization. In addition, the t ransfer correction shoul d 
be smaller as the spectrum is softer jAbel fc Haehneltll999D . 
Therefore, the background spectrum is unlikely to be as soft 
as a = 5 after the He II reionization. We should assume a 
harder index. 

A caveat against our results is in the correction factor 
of the transfer effect. Indeed, the effects of the photoelectric 
heating and the transfer correction can be offset. Thus, we 
should check this point by solving the cosmological radiative 
transfer in the future work. 

The current work may give a hint for the origin of the 
large b region of the observed IGM clouds. There is a possi- 
bility that the large b regions are localized. If the very small 
grains are located at the restricted areas, a systematic heat- 
ing owing to the dust photoelectrons may explain the large-6 
at the localized region. The principal role of the photoelec- 
tric heating in the IGM temperature also indicates that the 
theoretical equation of state in the IGM should be recon- 
structed by including the photoelectric heating effect. All of 
these works are being developed by the authors. 



ACKNOWLEDGMENTS 

We would appreciate valuable comments by the referee, Dr. 
A. Aguiree. We are grateful to R. Hirata for continuous en- 
couragement, and also thank H. Hirashita and S. Bianchi for 
very useful discussion. AKI is supported by the Research Fel- 
lowships of the Japan Society for the Promotion of Science 
for Young Scientists. 



REFERENCES 

Abel, T., & Haehnelt, M. G., 1999, ApJ, 520, L13 

Aguirre, A., 1999, ApJ, 525, 583 

Aguirre, A., & Haiman, Z., 2000, ApJ, 532, 28 

Becker, R. H., et al, 2001, AJ, 122, 2850 

Calzetti, D., 2001, PASP, 113, 1449 

Cen, R., 1992, ApJS, 78, 341 

Cheng, F. H., Gaskell, C. M., & Koratkar, A. P., 1991, ApJ, 
370, 487 

Cowie L. L, Songaila A., Kim T. S., & Hu E. M., 1995, AJ, 
109, 1522 

Crane, P., & Hoffman, A., 1973, ApJ, 186, 787 
Draine, B. T., & Lee, H. M., 1984, ApJ, 285, 89 
Draine, B. T., & Sutin, B., 1987, ApJ, 320, 803 
Eigenson, M. S., 1949, Soviet Astr., 26, 278 
Ferrara, A., Nath, B., Sethi, S. K, & Shchekinov, Y., 1999, 

MNRAS, 303, 301 
Hui, L., & Gnedin, N. Y., 1997, MNRAS, 292, 27 
Humason, M. L., Mayall, N. U., & Sandage, A. R., 1956, 

AJ, 61, 97 

Kitayama, T., Susa, H., Umemura, M., & Ikeuchi, S., 2001, 

MNRAS, 326, 1353 
Kogut, A., et al. 2003, ApJ, submitted ( |astro-ph/0302213F 
Laor, A., & Draine, B. T., 1993, ApJ, 402, 441 
Loeb, A., & Haiman, Z., 1997, ApJ, 490, 571 
Madau, P., Pozzetti, L., & Dickinson, M., 1998, ApJ, 498, 

106 

Nath, B. B., Sethi, S. K, & Shchekinov, Y., 1999, MNRAS, 
303, 1 

Nickerson, B. G., & Partridge, R. B., 1971, ApJ, 169, 203 
Perlmutter, S., et al., 1999, ApJ, 517, 565 
Riess, A. G., et al., 1998, AJ, 116, 1009 
Rowan- Robinson, M., Negroponte, J., & Silk, J., 1979, Na- 
ture, 281, 635 

Schaye, J., Theuns, T., Rauch, M., Efstathiou, G., & Sar- 
gent, W. L. W., 2000, MNRAS, 318, 817 
Sutherland, R. S., & Dopita, M. A., 1993, ApJS, 88, 253 
Spitzer, L., 1941, ApJ, 93, 369 

Spitzer, L., 1978, Physical Processes in the Interstellar 

Medium, Wiley, New York 
Takase, B. 1972, PASJ, 24, 295 

Telfer, R. C, Kriss, G. A., Zheng, W., Davidsen, A. F., & 

Tytler, D., 2002, ApJ, 579, 500 
Theuns, T., et al., 2002, ApJ, 574, Llll 
Todini, P., & Ferrara, A., 2001, MNRAS, 325, 726 
Weingartner, J. C, & Draine, B. T., 2001a, ApJS, 134, 263 

(WD01) 

Weingartner, J. C, & Draine, B. T., 2001b, ApJ, 548, 296 
Wright, E. L., 1981, ApJ, 250, 1 

Zheng, W., Kriss, G. A., Telfer, R. C, Grimes, J. P., & 
Davidsen, A. F. 1997, ApJ, 475, 469 



